Transverse momentum and centrality dependence of high-pr non-photonic electron 
suppression in Au+Au collisions at Ws NN = 200 GeV 
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The STAR collaboration at RHIC reports measurements of the inclusive yield of non-photonic 
electrons, which arise dominantly from semi-leptonic decays of heavy flavor mesons, over a broad 
range of transverse momenta (1.2 < pr < 10 GeV/c) in p+p, d+Au, and Au+Au collisions at 
^/s NN = 200 GeV. The non-photonic electron yield exhibits unexpectedly large suppression in central 
Au+Au collisions at high pr, suggesting substantial heavy quark energy loss at RHIC. The centrality 
and pr dependences of the suppression provide constraints on theoretical models of suppression. 
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High pt hadron production measurements at the Rel- 
ativistic Heavy Ion Collider (RHIC) show a strong sup- 
pression of the single-particle inclusive yields in nuclear 
collisions [3, [3, Q . The suppression is commonly thought 
to arise from partonic energy loss in dense matter due to 
induced gluon radiation with its magnitude depend- 
ing strongly on the color charge density of the medium. 
This makes it a sensitive probe of the matter created 
in heavy-ion collisions, where a quark-gluon plasma may 
form if sufficient energy density is achieved. 

Charm and bottom quarks are produced dominantly 
through high-Q 2 partonic interactions. Heavy flavor 
cross-sections and px spectra have been calculated at 
next-to-leading-order (NLO) for both p+p and A+A col- 
lisions [1, @, 0] , including nuclear matter effects Q ■ Al- 
though pQCD calculations agree well with heavy quark 
production in collider experiments at higher y /s [al , they 
disagree with recent RHIC measurements 0, flQj ]T Nev- 
ertheless, measurements of heavy quark production po- 
tentially provide new constraints on partonic energy loss 
mechanisms 



TABLE I: Corrections and systematic uncertainties for the 
non-photonic electron yield at pr = 2 and 8 GeV/c. 

Correction 



li, 12, i3, m \m m Uj 



Gluon radia- 
tion in a forward cone is suppressed for heavy quarks at 
moderate energy (dead cone effect) [III EH , with corre- 
sponding reduction in medium induced energy loss and 
less suppression of heavy-quark mesons than light quark 
mesons. 

Direct reconstruction of heavy flavor mesons via 
hadronic decay channels Q is difficult in the complex en- 
vironment of high energy nuclear collisions. Heavy quark 
production can also be studied through measurements of 
electrons (positrons) from semileptonic D and B decays. 
This Letter reports STAR Collaboration measurements 
of the non-photonic electron yield, (e + +e _ )/2, in p+p, 
d+Au, and Au+Au collisions at nucleon-nucleon center 
of mass energy ^fs^ = 200 GeV. The data extend sig- 
nificantly the px range of previous electron suppression 
studies [l8[ , to a region of phase-space where bottom de- 
cays are expected to be dominant. Large differences in 
energy loss are expected between c and b quarks in this 
region 14], and these measurements provide important 
new constraints on partonic energy loss mechanisms. 

STAR is a large acceptance apparatus comprising sev- 
eral detector subsystems within a 0.5 T solcnoidal mag- 
net field [l{|. The main detectors for this analysis arc 
the Time Projection Chamber (TPC) [20f and the bar- 
rel Electromagnetic Calorimeter (EMC)0. The EMC 
has a gas-filled Shower Maximum Detector (SMD) at a 
depth of ~ 5Xq to measure shower shape and position. 
A fast trigger based on single EMC tower energy enriches 
the electron sample at high pr- Electrons at moderate 
Pt were reconstructed from minimum bias and central- 
ity triggered Au+Au event samples, while EMC triggered 
events were used for px > 3 — 4 GeV/c. Au+Au data 
were divided into 3 centrality classes based on the track 
multiplicity measured at midrapidity. The integrated lu- 
minosity sampled by the EMC trigger is 100 nb" 1 for 



p+p 

2 GcV/c 8 GcV/c 



central Au+Au 
2 GcV/c 8 GeV/c 



Acceptance 0.84±0.05 0.75±0.15 

PID efficiency 0.25±0.03 0.50±0.03 0.13±0.03 0.45±0.03 

hadron contamination <0.01 0.20±0.04 0.03±0.03 0.22±0.05 

Bkgd. reco. cff. (en) 0.65±0.06 0.55±0.06 0.56±0.06 0.50±0.06 

Bremsstr. & Sp/p 0.86±0.14 1.05±0.05 0.9±0.1 l.lrtO.l 

EMC striggcr 1.00±0.08 1.00±0.05 

Cross section ±0.14 



p+p, 370 /ib _1 for d+Au and 26 /ib -1 for the most cen- 
tral Au+Au events. The charged particle acceptance is 
< T) < 0.7 and < <f> < 2ir, selected to minimize the 
radiation length of detector material interior to the EMC 
within the available EMC acceptance. 

The analysis has three main steps: selection of elec- 
trons; subtraction of background from decays and inter- 
actions in material; and residual corrections to the signal 
yield. Table U shows the major correction factors and 
uncertainties, which we now discuss in detail. 

Electron PID: Electron identification utilizes ionization 
energy loss (dE/dx) and track momentum from the TPC, 
together with energy and shower shapes from the EMC. 
Tracks with momentum p > 1.5 GeV/c are accepted if 
they originate from the primary vertex (distance of clos- 
est approach less than 1.5 cm) and project to an active 
EMC tower, with acceptance qemc ~ 75 — 85% of the 
EMC instrumented coverage. This reduced acceptance is 
due to dead or noisy electronics channels. Initial electron 
identification is based on p/E < 2, where p is the TPC 
track momentum and E is the energy of the EMC tower. 
Simulations show that this cut excludes ~ 7% of real 
electrons due to sharing of shower energy between tow- 
ers. Additional hadron rejection is based on the shower 
shape measured by the SMD. Figure QJi shows the dE/dx 
distribution for tracks passing the p/E and shower shape 
cuts. The curves show Gaussian functions fit to the dis- 
tribution, representing the yields of p + K, pions and 
electrons [22j. The parameters in the fit are the yields, 
widths, and overall dE/dx scale, with widths and the 
distances between centroids being quasi-free parameters, 
constrained by a model of energy deposition in the TPC 
gas [23|. 

Electrons are selected by cutting on TPC energy loss 
dE/dx min < dE/dx < 5.1 keV/cm. dE/dx min is around 
3.5 keV/cm, with the specific value having weak depen- 
dence on the event multiplicity and increasing slowly 
with track momentum, to optimize electron efficiency and 
hadron rejection while preserving more than 50% of the 
electrons in the dE/dx distribution. The residual hadron 
background satisfying the dE/dx cut is estimated based 
on Gaussian fits similar to those in Figure [TJ 

Table U shows the combined electron tracking and iden- 
tification efficiency ("PID efficiency"), determined by 
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FIG. 1: (a) dE/dx projections for 5 < pr(GeV /c) < 7 in cen- 
tral Au+Au events after EMC and SMD cuts. The lines are 
Gaussian fits for p + K, it, and electron yields, (b) Invariant 
e + e _ mass spectrum, (c) Ratio of inclusive and background 
electron yield vs. pt for p+p and Au+Au collisions. Vertical 
bars are statistical errors, boxes are systematic uncertainties. 



embedding simulated electrons into real events. It is sig- 
nificantly below unity due to tracking efficiency (~ 70%), 
exclusion of electrons due to the energy leakage to neigh- 
boring towers, and SMD response. Its increase from 
pr = 2 to 8 GeV/c is due to increasing SMD efficiency. 

Electron background: Background from photonic 
sources is due largely to photon conversions (~ 85%) in 
the detector material between the interaction point and 
the TPC (X/X ~ 4.5%) and tt° and <q Dalitz decays 
(~ 15%). The photonic electron yield is measured using 
the invariant mass distribution of track pairs detected in 
the TPC. One track of the pair is required to fall in the 
EMC acceptance, satisfying p > 1.5 GeV/c and electron 
PID cuts, with the other track having pt > 0.15 GeV/c 
within the TPC acceptance and a loose cut around the 
electron dE/dx band. Figure[T]D shows the invariant mass 
distribution of pairs with the same or opposite charge 
sign. The same-sign distribution is due to random (com- 
binatorial) pairs. An alternative combinatorial distribu- 
tion formed by embedding single simulated electrons into 
real events agrees with the same-sign distribution within 
statistical uncertainties. 

The shaded region in Figure [TJd is the difference be- 
tween the opposite and same-sign distributions and rep- 
resents the photonic yield. It exhibits a peak at zero 
invariant mass due to conversions, and a tail at non- 
zero mass due to Dalitz decays Selecting m < 150 
MeV/c 2 accepts ~ 98% of all tt° and r\ Dalitz pairs in this 
distribution. The efficiency Eb(pt) to identify a photonic 
electron in the EMC by this procedure was estimated by 
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FIG. 2: Non-photonic electron spectra. Vertical bars are sta- 
tistical errors, boxes are systematic uncertainties. The curves 
are scaled pQCD predictions for p+p Q. Cross section on 
right axis applies to p+p spectrum only. 



embedding [25| the main background sources (n° and 7) 
with a realistic momentum distribution derived from re- 
cent RHIC data (2|| into real events. 

The photonic electron yield N p h is calculated in each 
p T bin via N ph = (N unUke - N uke )/e B - Additional back- 
ground, mainly from u, <j>, and p decays, was estimated 
using PYTHIA [13 and HIJING [H simulations to be 
~ 2 — 4% of Nph 9j and is included in the systematic 
uncertainty of N p h ■ Figure [TJ; depicts the ratio of the 
inclusive to the photonic electron spectra for p+p and 
Au+Au collisions. The Figure shows a clear electron ex- 
cess. Within uncertainties, the non-photonic excess is 
independent of ccntrality at high pt- 

Non-photonic electron yield: The trigger efficiency was 
determined by comparing the electron candidate spec- 
trum in the minimum bias and triggered data sets. At 
high-pr the ratio of the spectra is compatible with the 
online scale-down factor applied to minimum bias events. 
The non-photonic spectrum is the difference of the in- 
clusive and photonic spectra. Additional corrections are 
applied for momentum resolution and brcmsstrahlung, 
determined from simulations. 

Systematic uncertainties: Systematic uncertainties 
were determined by varying cut parameters within rea- 
sonable limits. The uncorrelated systematic uncertainty 
of the electron yield is dominated by the electron iden- 
tification efficiency and photonic background reconstruc- 
tion at low pt and the correction for residual hadron 
background at high pt ■ 

Figure [2] shows the fully corrected non-photonic elec- 
tron spectra for 200 GeV p+p, rf+Au, and Au+Au col- 
lisions. The curves correspond to FONLL (Fixed Order 
Next-to-Leading Log) predictions Q for semileptonic D 
and B meson decays. The calculated spectrum is scaled 
by 5.5 (see below). 
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FIG. 3: Upper: ratio between measured non-photonic elec- 
tron yield and FONLL pQCD calculations Q for p+p colli- 
sions. Lower: relative contributions to FONLL distribution 
of c and b decays. 
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FIG. 4: The nuclear modification factor, Raa, for d+Au and 
Au+Au collisions at ^/s NN — 200 GeV. Error bars and uncer- 
tainties are described in text. 



Figure [31 upper part (points), shows the ratio of mea- 
sured to unsealed FONLL-calculatcd non-photonic elec- 
tron yield for p+p collisions. The calculation describes 
the shape of the measured spectra relatively well, though 
with a large difference in their overall scale. Better 
agreement is found at larger v / s_Sl- The same ra- 
tio is shown for published STAR Q and PHENIX 
measurements. The horizontal dashed line is at 5.5 ± 
0.8(stat) ± 1.7(sys), corresponding to the ratio between 
the total charm cross section measured by STAR Q to 
the central value predicted by FONLL 0, [I[ . The shaded 
band around that line shows the experimental uncer- 
tainty in this ratio. PHENIX data [10( exhibit a lower 
ratio and appear not to be consistent with the data re- 
ported here. The lower part (curves) shows the relative 
contribution to the FONLL calculation of charm and bot- 
tom decays, with the variation due to NLO uncertainties 
0, Li§| . The B-decay contribution is expected to be sig- 
nificant in the upper px range of this measurement. 

Modification of the inclusive particle production 
is measured by the nuclear modification factor [l[ 
(Raa(pt))- Raa is unity for hard processes without nu- 
clear effects. Figure 2] shows Raa{pt) for non-photonic 
electrons in d+Au and central Au+Au collisions. Error 
bars show the statistical uncertainties, boxes show un- 
correlated systematic uncertainties, and the filled band 
at unity is the overall normalization uncertainty. Raa 
for rf+Au is consistent with a moderate Cronin enhance- 
ment. Raa~ 0.2 for central Au+Au collisions at px > 3 
GeV/c, consistent with a previous measurement at lower 
Pt (1 81 ] - The suppression is similar to that for light 
hadrons at px > 6 GeV/c Q. 

Figure 2] shows predictions for electron Raa from scmi- 
lcptonic D- and B-mcson decay in central Au+Au col- 
lisions using calculations of heavy quark energy loss. 



Curve I uses DGLV radiative energy loss via few hard 
scatterings [LJ] with initial gluon density dN g /dy = 1000, 
consistent with light quark suppression. Curve II uses 
BDMPS radiative energy loss via multiple soft collisions 
with transport coefficient q. q is set to 14 GcV 2 /fm, 
though light quark hadron suppressi on p rovides only a 
loose constraint 4 < q < 14 GeV 2 /fm [15[. Both calcula- 
tions predict much less suppression than observed. 

This discrepancy may indicate significant collisional 
(elastic) energy loss for heavy quarks 13, 3(|. Curve III 
is a DGLV-based calculation including both radiative and 
collisional energy loss, together with path length fluctu- 
ations [l6[ . The calculated suppression is also markedly 
less than that observed. For Curve IV, the heavy quark 
energy loss is due to elastic scattering mediated by res- 
onance excitations (D and B) and LO t-channel gluon 
exchange 17( . This calculation also predicts significantly 
less suppression than observed. 

Dead cone reduction of energy loss is expected to be 
more significant for bottom than charm quarks in the 
reported px range. Curve V, which is the same calcu- 
lation as curve II but for D-mcson decays only, agrees 
better with the data. Since there is better agreement of 
data and theory for bottom than charm production at 
the Tevatron [8|, the scale factor 5.5 between calculated 
and measured p+p electron yields may overestimate the 
B decay contribution at RHIC, i.e. D decays may in fact 
dominate the electron yields in the reported px range, fa- 
voring calculation V. A direct measurement of D-mcsons 
at high-pr is required to understand energy loss of heavy 
quarks in detail. Finally, multi-body mechanisms may 
also contribute to heavy quark energy loss [31 1 . 

We have reported the measurement of high-px non- 
photonic electrons in p+p, d+Au, and Au+Au collisions 
at -%/s NN = 200 GeV. A pQCD calculation for heavy quark 
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production in p+p collisions undcrprcdicts the data, al- 
though it describes the overall shape of the pr distribu- 
tion relatively well. Large yield suppression is observed 
in central Au+Au collisions, consistent with substan- 
tial energy loss of heavy quarks in dense matter. The 
suppression is larger than that expected from radiative 
energy loss calculations, suggesting that other processes 
contribute significantly to heavy quark energy loss. This 
unique sensitivity to the energy loss mechanisms makes 
the measurement of heavy quark suppression an essential 
component of the study of dense matter. Full description 
of the interaction between partons and the medium will 
require further detailed measurements of charm and bot- 
tom separately. 
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